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ON THE DISTRIBUTION OF CUMULATIVE PARISIAN RUIN 


HELENE GUERIN AND JEAN-FRANgOIS RENAUD 


Abstract. We introduce the concept of cumulative Parisian ruin, which is based on the time 
spent in the red by the underlying surplus process. Our main result is an explicit representation 
for the distribution of the occupation time, over a finite-time horizon, for a compound Poisson 
process with drift and exponential claims. The Brownian ruin model is also studied in details. 
Einally, we analyze for a general framework the relationships between cumulative Parisian ruin 
and classical ruin, as well as with Parisian ruin based on exponential implementation delays. 


1. Introduction 

Historically, in classical actuarial ruin theory, the central topic has been the analysis of the 
probability of ruin, which is a measure to assess the overall risk of an insurance company. Ruin 
occurs when the surplus process falls below a certain threshold level for the hrst time, e.g. the 
insurer’s minimum capital requirement set by the regulatory body. Then, more sophisticated 
ruin-based risk measures using the timing and the severity of such a capital shortfall have been 
proposed, culminating with so-called Gerber-Shiu distributions; see mills]. 

However, as discussed in [9]: "[■■■] sometimes the event ruin has a very small probability and 
the portfolio is just one out of many existing in the company. The company can have enough 
funds available to support (or ask for external support) some negative surplus for some time, in 
the hope that the portfolio will recover in the future, allowing the company to keep this business 
alive. This can be regarded as an investment, since the process will recover in the future." This 
is why, in recent years, new ruin models, incorporating this realistic feature that a company is 
not immediately liquidated when it defaults, have been proposed. Therefore, new risk measures 
and models have been introduced: Parisian ruin (see, e.g., dllHllISllSo]), random observations 
(see, e.g., I2113I) and Omega models (see, e.g., 0, m and U)- In those papers, new dehnitions 
of ruin, bankruptcy and/or liquidation are proposed and studied. 

In this direction, we will introduce the concept of cumulative Parisian ruin, which is based on the 
time spent in the red by the underlying surplus process. While Parisian ruin has been inspired by 
Parisian options, the idea of cumulative Parisian ruin comes from cumulative Parisian options, 
another type of occupation-time options. In the definition of Parisian ruin (the duration of) 
periods of hnancial distress are considered separately, while cumulative Parisian ruin considers 
the aggregation of those durations. As expected, this notion of ruin is closely related to the 
occupation-time process of the underlying surplus process in the model. We will also analyze 
the relationships between cumulative Parisian ruin and classical ruin, as well as with Parisian 
ruin based on exponential implementation delays. 
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1.1. Insurance risk models. In this paper, we will mainly consider the following two classical 

actuarial rnin models: the Cramer-Lundberg model with exponential claims and the Brownian 
ruin model. First, in the Cramer-Lundberg model, the surplus of an insurance company/portfolio 
is modelled by a compound Poisson process with drift ,t> 0}, that is 

Nt 

(1) X^^ = X^^ + ct-^C„ 

i=l 

where c > 0 denotes the constant premium rate, N = {Nt,t > 0} is a Poisson process with 
rate A > 0 and the C^’s form a sequence of independent and identically distributed random 
variables (also independent of N) representing the claim amounts; each Ci is assumed to follow 
an exponential distribution with rate a > 0. In this model, the random variable Ci 

represents the aggregate claim payments made up to time t. 

On the other hand, in the Brownian ruin model, the surplus of the insurance company is modelled 
by a Brownian motion with drift > 0}, that is 

(2) Xf^ = X^^ + ct + aBt, 

where c > 0 denotes the constant premium rate, cr > 0 and B = {Bt,t > 0} is a standard 
Brownian motion. 

Note that both surplus processes X^^ and are special cases of a spectrally negative Levy 

process (SNLP), also known as a general Levy insurance risk process, that is a cadlag process 
X = {Xt,t > 0} with stationary and independent increments and no positive jumps. There is a 
rich literature on this family of stochastic processes for tackling classical ruin problems; see e.g. 

[mill]. 

1.2. Definitions of ruin and bankruptcy. Let X = {Xt,t > 0} be a surplus process. In 
what follows, the law of X with initial capital Xq = x is denoted by and the corresponding 
expectation by E^,. We write P and E when x = 0. 

The time of classical ruin is defined as 

Tg” = inf {t > 0: Xt < 0} , 

yielding the classical ruin-based risk measure P^; (rj" < oo), i.e. the probability of ruin. From 
a practical point of view, it is more interesting to consider the finite-time probability of ruin 
IPx (t'o" < t), where t > 0 is a fixed deterministic time corresponding to the time horizon of 
interest for the risk managers. 

In the last few years, in order to distinguish between default and bankruptcy, several concepts of 
ruin have been proposed. In particular, the idea of Parisian ruin has attracted a lot of attention; 
it considers the application of an implementation delay in the recognition of an insurer’s capital 
insufficiency. More precisely, it assumes that Parisian ruin occurs if the excursion below the 
critical threshold level 0 is longer than a pre-determined time called the implementation delay, 
or the clock. Two types of Parisian ruin have been considered: with deterministic delays or with 
stochastic delays. 

First, in IHIEQ], a Parisian ruin time (with a deterministic delay r > 0) is defined as 

Tr = inf {t > 0: t — gt > r} , 
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where gt = sup{0 < s < t : Xg > 0}, and studied for Levy insurance risk models. In other words, 
the company is said to be Parisian ruined the hrst time an excursion below zero lasts longer 
than the fixed implementation delay r. We will call this deterministic Parisian ruin, where 
deterministic refers to the implementation delay. 

On the other hand, in miiHiiis], exponentially distributed delays have been considered. A 
descriptive definition of Parisian ruin with exponential delay is given in [7] , for spectrally negative 
Levy processes. The idea is to mark the Poisson point process of excursions away from 0 by 
independent copies of a generic exponential random variable with rate g > 0. Then, the 
Parisian ruin time is defined as 

Kq = ini {t > 0: t - gt > , 

where gt = sup{0 < s <t : Xg > 0} and where each random variable is exponentially dis¬ 
tributed with rate > 0. We will call this exponential Parisian ruin. 

Note that for the two definitions of Parisian ruin mentioned above, there is a different clock for 
each excursion, either deterministic or stochastic. 

Remark 1.1. Note that the definition of exponential Parisian ruin corresponds to a specific 
choice of function oj{-) in an Omega model; see [^ 11211^ . 

The rest of the paper is organized as follows. In Section 2, after introducing the concept of 
cumulative Parisian ruin, we derive our main results. More precisely, in Section 2.1 and in 
Section 2.2, we obtain the distribution of cumulative Parisian ruin when the underlying surplus 
process is a Brownian motion risk process and a Cramer-Lundberg process with exponential 
claims, respectively. Section 2.3 is devoted to the study of the relationships between cumulative 
Parisian and classical ruin, as well as with exponential Parisian ruin. 


2. Cumulative Parisian ruin 


We now propose a new definition of actuarial ruin based on the occupation-time process of the 
surplus process X, namely the process 

t f l('_oo o)(^s)ds. 

Jo 

We dehne the time of cumulative Parisian ruin (at level r > 0) by 

(3) crr = inf|t>0: J l(_^o)(-^s)ds > J’j • 


Remark 2.1. Note that for cumulative Parisian ruin, the parameter r could be interpreted as the 
length of a clock started at the beginning of the first excursion, paused when the process returns 
above zero, and resumed at the beginning of the next excursion, and so on. In other words, it is 
the length of this single clock that is compared to the sum of all excursions below zero, as opposed 
to the previous definitions of Parisian ruin where a new clock is started at the beginning of each 
excursion. 


Remark 2.2. This definition of ruin is referred to as cumulative Parisian ruin due to its ties 
with cumulative Parisian options; see e.g. m Note also that this definition of bankruptcy has 
been used as a definition of default in structural credit risk models; see [22(125] . 
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We are interested in the following new ruin-based risk measure 

Fx (cTr < t) , 

where t > 0 is a fixed deterministic time. Clearly, for fixed values t,r > 0, we have 

> t} = < r| , 

so, the finite-time probability of cumulative Parisian ruin is given by 
Fx {ar < t) = 1 - Fx (^j l(_oo,o)(-^s)ds < 

(4) =1- Fx l(_oo,o)(-^s)ds G dv^ . 

Then, we also have 


(ct,. < oo) = 1 - Pa; l(_oo, 0 )(-’^s)ds < 

pr / poo \ 

J l(-oo,0)(^s)ds G duj . 


= 1 - 


As the distribution of cumulative Parisian ruin is closely related to the distribution of the 
occupation-time process associated with the underlying surplus process X over a finite-time 
horizon, we will now study the latter for the classical surplus processes. While there is a vast 
literature on the analysis of occupation times related to Brownian motion, it turns out that little 
is known about the distribution of the occupation-time process for a compound Poisson process 
with drift over a finite-time horizon. 

Of course, it is very difficult to obtain analytical expressions for the probability of finite-time 
cumulative Parisian ruin for a general insurance risk process. However, there is a quite simple 
expression of its Laplace transform when the initial capital is zero in the case of spectrally 
negative Levy processes. 


Proposition 2.1. If X is a SNLP such that E[Aii] = ip'{0+) > 0, then, for p,q > 0, 


( 5 ) 


f 


-P^E 


^- 9/0 l{-oo, 0 ) (W)ds 


dt 


^{p + q) 1 
p + q $(p)’ 


where ^{q) is the largest root of Lundberg’s equation — q = 0, for a fixed q > 0, and where 
-iP{6) :=ln (E [e^^i]). 


Proof. See e.g. Remark 4.1 in |19) . 


□ 


Since 


f 


■P*E 




dt = 


(^'l(-oo,0)(^.)ds G d^r^ dt. 


we obtain the distribution of the occupation time by inverting this double Laplace transform, 
which is our objective for the two main ruin models considered in this paper. 
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2.1. Brownian ruin model. One of Paul Levy’s arcsine laws gives the distribution of the 
occupation time of the positive/negative half-line for a standard Brownian motion B = {Bt,t > 
0}. More precisely, for 0 < s < t, 


P 



l(_oo,0)(-B«)d« < s 



This result was then extended to a Brownian motion with drift by Akahori [T] and Takacs |24| . 
Recall that in the Brownian ruin model of Equation ([2|), the surplus is given by 

+ ct + aBt, 


1 


-|- 2a^q — c 


where c > 0 and cr > 0. In this case, = c9+^6^, with right-inverse ^(g) = 

From Proposition 12.11 the double Laplace transform of the cumulative Parisian rum cr,- is given 
by, for p,q>0, 





dt 


y/c? X 2a‘^{p -h q) - c 1 

P + Q (? X 2a^p — c 


We recognize (1.4) in [T] (given for the case a = 1) and then, by Laplace inversion, we recover the 
generalized arcsine law of the occupation time of the negative half-line for a Brownian motion 
with drift. 


Theorem 1 (Akahori (1995) [I]). For a fixed t > 0, the distribution of the occupation time of 
the negative half-line, when X^^ = 0, is given by 


P 


l(-oo,0)(^«'^)du G ds ) = 


2 ae 


-(c2/2(t2)s 


CT^ 


V2 


- cN 


ITS 


c^/s 


X <0 + 


ae 


a 

-{cy2a^){t-s) 


y/2Tr{t - s) 


-cN 


Zy/t — 


a 


ds. 


for d < s <t, where N denotes the tail of the standard normal distribution, that is 




dy. 


Clearly, using Equation ([3]), we have an explicit expression for P (<7^ <t) in the Brownian 
model. Using the continuity of Brownian motion, we can easily obtain an explicit expression 
for Pa; (cJr < t), i.e. when = X > 0. Indeed, using the strong Markov property and the 

continuity of we easily deduce that, for s < t, 


P. (^*l(_oo,o)(^™)dn<s) 

/ t — S / pt — U \ 

Fx (tq- G dn) P yj l(_oo, 0 )(-^™)du < sj du, 
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where it is well-known (see e.g. |15) ) that 


Pcc (tq G du) 


V27ra^u^ { 


(x — cu)^ 
2a‘^u 


du. 


From the above, we easily deduce the following expression for the probability of cumulative 
Parisian ruin: 


Pa; (fJr < t) = 1 




(x — cu)‘^ 
2a‘^u 


du. 


2.2. Cramer-Lundberg model with exponential claims. Recall that in the Cramer-Lundberg 
model with exponential claims introduced in Equation ([T]), the surplus is given by 

Nt 

i=l 

where c>0,A^ = {A^t,t>0}isa Poisson process with rate A > 0 and the Cj’s form a sequence of 
independent and identically distributed random variables following an exponential distribution 
with rate a > 0. 

In this model, for p > 0, Lundberg’s equation takes the form 'ijj{6) —p = 0, where 
and its solutions are given by 

^(p) = ^ ^ + V^) > 

0(p) = ^ (^p + A - CO - v^) j 

where Ap = (p + A — ca)^ + 4cap. Note that we can also write Ap = (p + A + ca)^ — 4caA. 
Clearly, for a fixed t > 0, the distribution 

has a mass at at 0, which is the survival probability P (r^ > t). We now derive from Proposi¬ 
tion 12.11 an expression for this distribution. 


Proposition 2.2. For a fixed t > 0, we have 


( 6 ) 

with 

( 7 ) 


P 


(^J e ds^ = at6o (ds) -h ats (A - ca(l - a*)) l(o,t)(s)ds. 


at= [1- 


ca 


TT 


/: 


VT^ 


-1 A -|- CO -|- 2 \/caX 




u 


and z_|_ := max(z, 0). 


Note that there exists other equivalent expressions for a* = P (tq > tj in the Cramer-Lundberg 
model with exponential claims. See e.g. |in) . or [H Proposition 1.3] and the references therein. 
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Proof. For simplicity, we will use X instead of in the proof. From Proposition 12.11 we know 
that the double the Laplace transform of the occupation time of the negative half-line is equal 
to 


( 8 ) 



Q-lfo l(-oo.O)(^s)ds 


dt 


^{p + q) 1 

p + q 


Since ^{v)/v converges to 1/c when v goes to oo, we deduce that ([8]) converges to (c<k(p)) ^ 
when q goes to oo. As mentioned above, this is due to the fact that has a 

mass at 0, which means that its distribution under P can be written in the following way: 

P l(_oo,o)(^n)d7t G ds^ = at5o(ds) ht(ds), 


where at G [0,1] and ht{ds) is a nonnegative measure on [0,t] with no mass at zero. 
We easily see that 


( 9 ) 

Then, by linearity. 


poo poo 

Ca{p) := / e~^^atdt = / e“^*P {tq > t) dt = 

Jo Jo 

the double Laplace transform of {t, s) i— ht{ds) must be equal to 


c$(p)' 


( 10 ) 


Ch{p, q) 


1 / c^{p + q) _ 

c^(p)V P + q ) 

1 / c9{p + q) \ — ca\ 

c$(p) \ p + q P + q ) 
1 / a A — ca \ 

c<^{p)\^{p + q) p + q)' 


where the second equality is obtained using the explicit expressions for <k(-) and 6{-), while the 
third equality is a consequence of the following identity: c<k(p -|- q)0(p + q) = —ot{p + q). 

Using basic properties of Laplace transforms and Equation ([9]), the expression of the Laplace 
transform of (t, s) e-)• /if(ds) implies that 


ht{d.s) = at-s{caas + A - ca)l(o,t) (s)ds. 

Consequently, to identify the law of l(^_^ Qj(Xu)du, we need to identify at- Instead of using 
one of the available expressions from the literature, we derive a new expression. 

First, using the conjugate expression of square roots and the explicit form of <k(p), we notice 
that 


Ca{p) = 


y/{p + X + ca)^ — dcaX — (p -|- A -|- ca) 


2cap 

-2A 


1 

P 


P 


p + \ + ca)‘^ — 4cq:A -|- (p + A -|- ca 


+ 


P 


As, for functions / and g, Cf{p) = ^ Jq J^q(q)dq implies f(t) = s ^g{s)ds and the Laplace 
transform of t i— )• Ii{kt) is equal to p i— )• k(^\/p^ — kJ + {jp^ — k‘^) with Ii the modified 
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Bessel function of the first kind of order 1 (see e.g. |23t Supplement 4] for more details), we 
deduce that, when A > ca, or equivalently E[Xi] < 0, 



and, when A < ca, or equivalently E[Xi] > 0, 


(^+“)*/i(2V^s)ds, 


A / \ /*oo 1 

at = l - + \ — / (2 s)ds. 

ca \ ca s 

At last, the final result is obtained using Fubini’s theorem and the following integral form for 
the modified Bessel function: for s > 0, 


h{s) 


-f 

TT J-1 


e-"“Vl - uHu. 


□ 


As opposed to the Brownian case, it is more delicate to obtain the corresponding distribution 
under P^, for x > 0, i.e. when Xq^ = x > 0, since X^- < 0 a.s.. 

Again note that, for a fixed t > 0, the distribution 

has a mass af at 0, which is the survival probability P^, (tq~ > t). 


Theorem 2. For a fixed t > 0, we have 

rt 


(11) Px l(_oo, 0 )(^“)d^x G ds^ = of(5o(ds) + (of_^ + kf_^) (A - ca (l - a°)) l(o,t)(s)ds. 


with 


= 1 - Ae"""^ / 


cay/s{s + x/c)^ 


ds. 


( 12 ) 


kf = e-"^ - 1 


+ xaXe~°‘^ J Iq ^2-\/caAY^s(s + x/c)'j — /2 ^2\/caAv^s(s + x/c)^ 


ds, 


where represents the modified Bessel function of the first kind of order v, which can be written 
as: for s > 0, 


Io{s) = 



and l 2 {s) 


e-"“(l-u2)^/^du. 
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Proof. For simplicity, we will use X instead of in the proof. First, note that we have the 
following relationships between <h(-) and for p, g > 0, 

(13) $(p) - 6»(p) = 

(14) ci,(p)0(p) = _^, 

c 

(15) {0{p + q) - ^{p)){0{p + q) - 0{p)) = -^{9{p + q) + a). 

From m Corollary 1], for x G M and p,q > 0, we deduce an explicit expression for the double 
Laplace transform of the occupation time of the negative half-line 



Q-lfo l(-oo.0)(^dds 


dt 



^{p + q) - ^(p) 


^(P+q,-g)(^)^(p,g)(_y) _ >v(P>'?)(x - y) 


dy, 


where, in the case of a Cramer-Lundberg risk process with exponential claims (see e.g. |13| for 
more details), for x > 0, 

y(p^<i)(x) = g \ a + Hp + q) _ _^±0ip±^ e{p+q)x 

y^g[Hp + q)-Hp) e{p + q)-Hp) \ 

and 

y{p+q-q)(j.) = g ( «+ ^(P) $(p)x _ «+ Gjp) 9(p)x\ 

V^(p + g) - ^(p) ^{p + q) - 0{p) / 

and where, for x < 0, = exp(<l>(p)x) and = exp (4)(p-|-( 7 )x). Also, for 

X > a. 


>V^P’«)(x) 


^ a + ^{p + q) 

\/ ^p+q^p 

_ ^a+_0(p_-^^0(p+g)(,j,_„) 

\/ ^p+q^p 


a + <I>(p) 

^{p + q)-^{pf 

a + <l)(p) 

^(p + g) - 4>(p) 


$(p)a 


,g'*’(p)a 


« + ^(P) ^Biv^a 

^{p^q)-0{p) 

_ « + ^(P) e(p)a 

(9(p + g) - 6»(p) 


for X G [0,a], >Va^’'^^(x) = W^\a), and finally, for x < 0, W^'^\x) = 0. Therefore, for p < 0, 
we have 


W^J>’^\x-y) 


a + <l>(p -b g) 


\J ^p+q 


y{p+q,-q) 


^ a + 9{p + q) [■ _a+_$(p)_^,j,(p),j, 

^ i/Ap+qAp |_6'(p + g) - ^(p)^ 


a+ 6{p) e{p)x 

<9(p + g) - 6'(p) 


Q-S{p+q)y _ 
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We divide the computation of the integral into three parts. We first consider the integral on the 
negative half-line. For y < 0, using ffT5|l and (I13p . we have 




a -I- 6{p + q) 

\/ ^p+q 


^(P + g) ^{p) y(pJ^q-q), 

e{p + q)-m 


+- 


a + ^(p) ^Mp)x _ a + 0{p) „e(p), 


v^V^(p + g) - ^(p) 


6{p + q) - 9{p) 


c{a -h 6{p)) 

L 

_ a + 9{p) 
c{^{p + q) - 9{p)) 

Consequently, 

/■O 

( 16 ) / 


(Cl>(p + q)- - {9{p + q)- m) 


^(p + g) - 9{p) 


^-e{p+q)y 


^e{p)x^-e{p+q)y 


Qd{p)x^-e(p+q)y ^ 


£(P±9)-_ y) 

g J 


dy 


a d- 9{p) 


J{p)x 


c9{p + q){^{p + q) - 9{p)) 

We now compute the integral on the positive half-line for each term. We have 

(17) ^ ^(p+g)-^(p) ^ ^{p+q,-q) (_y)dy 

g 7o g^(p) 

At last, we notice that x — y<x when y > 0 and then —y) = —y). Consequently, 

(18) - r wiP^'i\y)dy = - rw(P)(y)dy = - ^ . . . 

Jo Jo Hp)y^ ^(P)V^ c^p)9{p) 

At last, summing up the three integrals in (1161) . (1171) and (I18D . we obtain 


(19) 



e“'?/o ^(-°o,o)(W)ds 


dt = 


{a + 9{p))^{p + q) 
c($(p + y) - 9{p)) 


1 

$(p)6»(p) 


1 

^>(p + g)6»(p + y) 




a 

c$(p)6»(p)' 


The rest of the proof is similar to the proof of Proposition 12.21 i.e. we want to invert the double 
Laplace transform in (HU). Fix p > 0. Since when q ^ oo we have $(p -|- g) —>• oo and 
9(p + q) ^ —a, then the previous quantity converges to 


( 20 ) 


a+ 9{p) ^ a + 9{p) ^ 1 

c^{p)9{p) c^{p)9{p) ap p 
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Since the distribution of has a mass at 0, 


P, 


l(_oo, 0 )(^n)du G ds^ = of(5o(ds) + hf(ds), 


where af = Pa;(T'cr > i) S [0)1] hf(ds) is a nonnegative measure on [0,t] with no mass at 
zero. From (j20l) . we deduce that the Laplace transform of t af is equal to 1 — q^{p)^ 

and then, from (fT^ . the Laplace transform of {s,t) i—)• /if(ds) is 


CH^{p,q) = 


{a + e{p))^{p + q) 
c{^{p + q) - e{p)) 


J)(p)x _ Q + 0{p) ^B{p)a 

c^{p)6{p) 


/^{p)e{p) ^{p + q)e{p + q) 

Taking x = 0 in (1191) and comparing with ([S]), we deduce that the Laplace transform of {t, s) 
hf{ds) is equal to 

qS(p)x 


( 21 ) 


Ch^{p,q) = 


c$(p + q) 

p + q 


- 1 


c$(p)' 


Therefore, to identify the law of Jq l(_oo^o)(-^«)dii, we need to identify 

at = Pj, (to" >t) = 1- [ Pj, (to" G ds) . 

Jo 

Using the explicit expression for P^, (r^ G ds) given in cni Equation (D.17)], we deduce that 


a? = 1 - Ae- 


f 


^-{\+ca)s 


fo ^2\/AmA/s(J+”x7^^--j-—^/2 ^2\/Aca-\/s(J+"x7^^ 


ds. 


where represents the modified Bessel function of the first kind of order zz, which can be written 
as 

,21-1/2 


Io{s) = — I e ™(l —tt"') "'"du 


3-^(1 _ ^2^ and his) = ^ r e-'*“(l - u^^^^du. 

J~1 uTT 

Moreover, using ()13p . the Laplace transform of can be written in terms of the Laplace trans¬ 
form of a^: 

c^ip + q) 


Ch^{p,q) = 


p + q 


- 1 




We then deduce that (from the proof of Proposition 12.21) 

/if (ds) = ca ^ - 1^ (af_^ -b kf_s) l(o,i)('S)ds, 

with Ck^ip) = — l). As, for functions / and g, Cfip) = 1 Cgiq)dq implies /(/) = 

/o s“^ 5 '(s)ds, we have kf = s“^mfds -|-e“"® — 1, with Cm+{p) = —x6'ip)e^^^^. Let us now 
explicit the function mf. We have 

—ca—\J (p-t-A-t-ca)^—4 c«a) 

— (p+-^+cap—4caA—(p+\+ca)') 


£m^(p) = 1 - 


p -b A -b ca 


xe 


2c 


(p -b A -b ca)^ — 4caA 
p -b A -b ca)2 — 4caA — (p -b A -b ca)^ 


y^(p -b A -b ca)2 — 4caA 
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Using similar inversion techniques as in [lOl (D.13)-(D.15)], the inverse Laplace transform of the 
previous expression can be expressed in terms of the modified Bessel function Ii : 

j, (2V^v' t(t + ^/c) ) ■ 


= 


2c ^yt{t + x/c) 

Since Ii{z) = z{Io{z) — we have 


mf = xaAe-"^“e-(^+'^“)* 


Io(2Vcc^^/t{^+lcJ^'J — I 2 ^2\/ ca\-\/t{t + x/c)J 


and then 


kf = e-“^ - 1 


+ xaAe 


/ 


g-l^-(X+ca)s 


Iq ^2\/ caX^ys{s + x/c)^ — I 2 ( 2 VcaXy^s(s + x/c)^ 


ds. 


□ 


If we set a; = 0 in the definition of kf in m, then clearly = 0. Therefore, Theorem [2] is a 
generalization of Proposition 12.21 


2.3. Relationships with other definitions of ruin. When r —)■ 0, we should intuitively 
recover the classical probability of finite-time ruin. Indeed, when r is very small, then almost no 
grace period is granted when the surplus goes below zero. 

First, for a: > 0 and t > 0, we have for every r > 0 that 

rt 


l(_oo, 0 )(^ 6 )ds > r| C {tq < t} 


so 


Fx (Cr < t) < Fx (tq < t) . 

Therefore, as a risk measure, the probability of cumulative Parisian ruin is less conservative than 
the probability of ruin. 

We can obtain a much stronger relationship between cumulative Parisian ruin and classical ruin 
in the full generality of a Levy insurance risk model. Recall that a Levy insurance risk process is 
a process X = {Xt,t > 0} with stationary and independent increments and no positive jumps. 

Proposition 2.3. Assume X is a Levy insurance risk process. For all x > 0, the time of 
cumulative Parisian ruin ar eonverges ¥x-almost surely, as r —)• 0, to the time of elassieal ruin 

^0- 

Proof. First, we show convergence in probability (with respect to F^,, for an arbitrary value of 
X > 0), that is, for every e > 0, 

Fx(|a^-ro-| >e) ;:^0. 

Indeed, for any e > 0, we note that ar > Tq and, by the strong Markov property. 


Fx(|o-r - Tq I > e) = Pa; 


re+Tn 


l(-oo,0)(-’^s)ds < C = Ea 


Fx^.ry l(_oo,0)(^5)ds < c 
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Therefore we want to show that Px _ (Jq o)(^s)ds = O) 

L '^0 

since the process is skip-free upward, 

Px^_^y l(_oo,o)(-^^f)ds = 0^ = 0. 


0. Clearly, on {X^- < 0}, 


On the other hand, if x > 0, it is well known that P^, = Oj if and only if X has a Brownian 

component. Finally, recall that a SNLP has paths of unbounded variation (e.g. with a Brownian 
component) if and only if 0 is regular for (—oo,0), which means that when the process sits at 
zero it will enter (—oo,0) immediately (and will start to accumulate occupation time). Putting 
the pieces together concludes the proof of convergence in probability. 

Since r i—)• cJr is a non-decreasing function, we deduce that cr^ converges to Tq a.s. when r goes 
to 0. □ 


At first sight, it is very tempting to think of exponential Parisian ruin as being one Laplace 
transform away from deterministic Parisian ruin. As discussed in |2n| . due to the choice of 
different exponential clocks for different excursions below zero, this is not the case. In fact, it 
turns out that cumulative Parisian ruin is the notion of Parisian ruin that can be connected to 
exponential Parisian ruin through a Laplace transform, as the next result will show. 

Note that the next result holds for any Levy insurance risk process, including of course the 
Cramer-Lundberg risk process and the Brownian risk process studied previously. 


Proposition 2.4. Assume X is a Levy insurance risk process and let be an independent 
exponential random variable with rate q > 0. Then, for all initial surplus x > 0, Hq and (Te^ are 
such that, for all x > 0, t > 0, 

Ex (Kg < t) = Ex < t) 

and, moreover, 

Ej; {Kq < oo) = Ea; < Oo) . 


Proof. Comparing [T3l Corollary 1] and |71 Corollary 1.1], we observe that, for all y E (—oo,0], 

r t 1 r 

q / e-®*Ex dt = E^ e-'^'L X,, £ dy, Kg < oo 

Jo ^ ^ ^ 

Using Fubini’s theorem and then integrating over (— 00 ,0] (with respect to y) both distributions, 
we get 


/•OO 

/ e-®*E, 

Q-‘!fo'^(-oo,o){Xs)ds.x^ < 0 

dt = Ex 

e < 00 

Jo 



- 


Note that we can write 


Ex 




, Kg < OO 


f 


Therefore, by Laplace inversion, we have, on (0,oo), 


gEx 


3-'?/o l(-oo,0)(^s)ds. ^ g 


(Kg G dt) . 

dt = P3, (Kq £ df) . 
















14 


Integrating over (0, t) for an arbitrary value of t yields 


P, 


{Kq<t)= [ qE^ 

Jo 


ds 


= Ex 

= Ex 

= 1 - E, 


[ ?! (-00,0) )e"'' foh-^.o) 

Jo 

I _ e“'?io i(-oo,0)(^s)ds 
e^'j/o l(-oo,0)(^dds 


So, we have 

(22) Px (Kq < t) = 1 - Px l(_oo,0)(-’^s)ds < = Pa; ((Te,j < t) , 

and, in particular, 

(23) Px (/tg < oo) = 1 - Ea 


,-q/o°“l(_oo,0)(^dds 


= Ex (<Te,j < oo) 


□ 


Remark 2.3. The previous proof highlights a strong relationship between the distribution of the 
occupation-time process and exponential Parisian ruin, namely identity (|22P The equality in 
Equation (1231) was first discussed in |18| ; note that Equation (I22p is a finite-time version of it. 
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